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ABSTRACT 

We have obtained FUV spectra of two clusters of galaxies with FUSE, the Far Ultraviolet Spectroscopic 
Explorer. The Coma cluster was observed for a total of 28.6 ksec, the Virgo cluster for 10.9 ksec. Neither 
spectrum shows significant O VI AA1032, 1038 emission at the cluster redshift. Such emission would be 
expected from the warm [(5 — 10) x 10^ K] component of the intracluster medium (ICM) that has been 
proposed to explain the excess EUV and SXR flux present in EUVE and i? 05*^4 T" observations of these 
clusters. Our 2-a upper limits on the O VI A1032 flux from Coma and Virgo exclude all published 
warm-gas models of the EUV excess in these clusters. 

Subject headings: galaxies: clusters: general — ultraviolet: galaxies 



ter place tight constraints on the emission integral of any 
warm component of their ICM. 

2. OBSERVATIONS AND DATA REDUCTION 

FC/Si? comprises four separate optical systems. Two em- 
ploy LiF optical coatings and are sensitive to wavelengths 
from 990 to 1187 A, while the other two use SiC coatings, 
which provide reflectivity to wavelengths as short as 905 A. 
The four channels overlap between 990 and 1070 A. For a 
complete description of FUSE, see Moos et al. (2000) and 
Sahnow et al. (2000). 

The FUSE spectrum of the Coma cluster was obtained 
in 17 separate exposures on 2000 June 18 and 19. Each 
was centered on 12'^59'"49!0, +27°57'46" (J2000, or I = 
57.61,6 = H-87.96 in Galactic coordinates), near the clus- 
ter center. The total exposure was 28608 s, with 23553 s 
obtained during orbital night. We use the entire 29-ksec 
data set in our analysis. Our spectrum of the Virgo clus- 
ter is a combination of data from two locations near the 
cluster center. Two exposures, centered on 12'^31™07;3, 
-fl2°23'46" (/ = 284.03,6 = -1-74.52) and totaling 2242 
s, were obtained on 2000 June 13. Seven exposures, cen- 
tered on 12i^31'"13M, -|-12°22'10" {I = 284.2,6 = 74.50) 
and totaling 8688 s, were obtained on 2000 June 17. The 
total integration time is 10,930 s, all of it during orbital 
night. All observations were made through the 30" x 30" 
(LWRS) aperture. 

The standard data-reduction software are used to 

confirm that none of our data were obtained during pas- 
sages through the South Atlantic Anomaly or at low earth- 
limb angles, to separate day and night observations, and 
to confirm that the dead-time correction for both obser- 
vations is negligible. To reduce the detector background, 

^ Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE is operated for NASA by the Johns 
Hopkins University under NASA contract NAS5-32985. 

Current address: Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA 
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1. INTRODUCTION 

Recent observations with the Extreme Ultraviolet Ex- 
plorer (EUVE) provide evidence that a number of clus- 
ters of galaxies are strong emitters of EUV radiation (Lieu 
et al. 1996a,b; Bowyer et al. 1997; Mittaz et al. 1998). This 
emission is substantially brighter than would be expected 
from the well-known. X-ray-emitting intracluster medium 
(ICM; Forman & Jones 1982). Signatures of this "soft ex- 
cess" are sometimes present in the lowest energy-resolution 
band of the i?05'^T' PSPC, but its presence remains con- 
troversial (Arabadjis & Bregman 1999; Drake 1999; Krick 
et al. 2000). At present, two groups agree that an EUV 
excess exists in both the Virgo and Coma clusters, but 
they differ on the intensity and morphology of the emit- 
ting regions (Lieu et al. 1999; EnBlin et al. 1999; Bowyer & 
Berghofer 1998; Bowyer et al. 1999; Berghofer et al. 2000). 

The EUV excess in clusters of galaxies was originally at- 
tributed to a diffuse, (5 — 10) x 10^ K thermal gas compo- 
nent of the ICM. Gas at these temperatures cools rapidly, 
however, requiring cither a substantial mass of cooling gas 
or a heretofore unknown energy source to reheat it (Fabian 
1996). If a reservoir of warm gas were present in the cores 
of clusters of galaxies, it would emit strongly in the far- 
UV (FUV) resonance lines of O VI AA1032, 1038 and C IV 
AA1548, 1551 as it cools through temperatures of a few 
times 10^ K (Edgar & Chevalier 1986; Voit et al. 1994), 
yet Dixon et al. (1996) were unable to detect either dou- 
blet in the spectra of five clusters of galaxies obtained with 
the Hopkins Ultraviolet Telescope (HUT). To probe more 
deeply for this emission, we have observed the Virgo and 
Coma clusters with the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE). Our limits on the O VI flux of each clus- 
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photon events with pulse heights less than 4 or greater 
than 15 (in standard arbitrary units) are excluded from 
the data set. Subsequent data reduction is performed by 
hand. Photon events from regions of the detector with 
known defects arc removed. The standard extraction win- 
dow defined for the LWRS aperture is centered (in the 
dimension perpendicular to the dispersion axis) on the dif- 
fuse airglow emission lines. 

While the resolution of the FUSE spectrograph is ap- 
proximately 15 km s,'^ for a point source, diffuse emission 
filling the LWRS aperture yields a line profile that is well 
approximated by a top- hat function with a width of ~ 106 
km s~^. Our spectra are quite faint, so we bin the Coma 
data by 8 detector pixels, or about 15.5 km s^^. We bin 
the data from the shorter Virgo observation by 16 pixels. 
We distinguish between detector pixels and 8- or 16-pixel 
bins throughout this paper. 

We omit the final steps of the standard FUSE data- 
reduction pipeline for several reasons: first, the pipeline's 
corrections for differential Doppler shifts and grating and 
electronics drifts total less than 10 detector pixels in the 
dispersion direction. Second, the distortion correction, 
applied independently for each of the instrument aper- 
tures, can cause background photons to "pile up" in the 
region between two nominal extraction windows, forming 
a horizontal stripe in the processed detector image. These 
stripes occasionally overlap our preferred extraction win- 
dow. For bright point sources, this effect is unimportant, 
but for faint, diffuse sources, we prefer not to risk including 
such features in our extracted spectrum. Finally, we wish 
to include the detector background as a free parameter in 
our spectral models, so forgo the background subtraction 
steps included in the standard pipeline. 

The FUSE flux calibration, based on theoretical models 
of white-dwarf stellar atmospheres, is believed accurate to 
about 10% (Sahnow et al. 2000). Corrections to the nomi- 
nal FUSE wsNelengih scale are derived from the measured 
positions of airglow features in each detector segment and 
are good to about 0.01 A. In practice, only the portion of 
the Coma spectrum obtained during orbital day contains 
airglow features strong enough to allow an independent 
wavelength correction. We thus apply the Coma correc- 
tion to the Virgo data; the resulting wavelength solution 
is consistent (within the errors) with the positions of mea- 
sureable airglow lines. Error bars are assigned to the data 
assuming Gaussian statistics, then smoothed by 9 bins to 
remove small-scale features in the error spectrum with- 
out significantly changing its shape. Segments of the flux- 
and wavelength-calibrated LiF lA spectra, showing O VI 
AA1032, 1038 at the redshift of each cluster, are presented 
in Fig. 1. We assume z = 0.0231 for Coma (Struble & 
Rood 1999) and z = 0.0036 for Virgo (Ebeling et al. 1998). 
Diffuse emission from our own Galaxy is present in both 
spectra; these results are presented in Dixon et al. (2001). 

3. SPECTRAL ANALYSIS 

In our analysis, model spectra are fit to the flux- 
calibrated data using the nonlinear curve-fitting program 
SPECFIT (Kriss 1994), which runs in the IRAF^ environ- 
ment, to perform a minimization of the model param- 



eters. We find that, because of the wide disparity in the 
effective area of the various detector segments, combin- 
ing spectra from different segments does not significantly 
improve the signal-to-noise ratio. We thus use data only 
from the segment with the highest effective area at the 
wavelength of interest. 

To set limits on the flux of an emission feature, we first 
fit a linear continuum to the spectrum. We then add a syn- 
thetic emission feature to the data array, raising its flux 
until for the best-fit model with an emission feature 
differs from for the best-fit model without an emission 
feature by Ax^ = 9 (corresponding to a 3-a deviation for 
one interesting parameter, in this case the flux in the emis- 
sion line; Avni 1976). We quote as a S-cr upper limit the 
flux of the model emission line that best fits the synthetic 
emission feature originally added to the data. Upper limits 
derived in this way are unaffected by small-scale features 
in the observed continua. The best-fit linear continuum 
and our 3-a upper limits to the flux of the O VI doublet 
for each cluster are shown in Fig. 1. 

The observed profile of a diffuse emission feature repre- 
sents a convolution of its intrinsic profile with the 106 km 
s~^ top- hat function discussed above. In our models, we 
assume the intrinsic profile to be a Gaussian with FWHM 
= 40 km s~^, corresponding to the thermal width of a 
5 X 10^ K gas. Broader lines, perhaps due to bulk motions 
in the gas, would yield higher upper limits, as the flux is 
spread over a greater region of the detector. For c;xaniple, 
lines with an intrinsic FWHM = 200 km s~^ yield upper 
limits approximately 70% higher than those of 40 km s~^ 
lines. 

No significant FUV line emission is seen at the cluster 
redshift in either spectrum. Our limits on the fluxes of 
various lines predicted to be strong in thermal-gas (Lan- 
dini & Monsignori Fossi 1990), mixing-layer (Slavin et al. 
1993), and shock (Hartigan et al. 1987) models are pre- 
sented in Table 1. Dereddened intensities are calculated 
assuming the extinction parameterization of Cardelli et al. 
(1989, hereafter CCM) with E{B-V) = 0.008 for Coma 
and 0.030 for Virgo (Schlegel et al. 1998) and Ry = 3.1. 
Uncertainties in the reddening correction are discussed in 
Dixon et al. (1996). 

At the wavelength of redshifted O VI in the Coma clus- 
ter (1055.77, 1061.59 A), the LiF lA detector channel has 
an effective area Aes = 24.3 cm~^, by far the highest of 
the four FUSE channels (Sahnow et al. 2000). Unfortu- 
nately, a bright band of scattered light contaminates the 
LiF lA spectrum from about 1045 to 1058 A. Present even 
at night, the band is clearly visible in Fig. 1. It shows con- 
siderable structure, including an apparent emission feature 
at 1054.5 A. Because the intensity of the band varies on 
small spatial scales, background spectra extracted from de- 
tector regions above and below our spectral window cannot 
be used to predict its shape in our spectrum. The statis- 
tical significance of the 1054.5 A feature is about 3.3 cr; if 
real, it should be present in the spectrum from the less- 
sensitive {Af.g = 15.5 cm~^) LiF 2B channel at a signifi- 
cance of > 2 CT. No such feature is seen, and we conclude 
that it is an artifact of the stray-light stripe in the LiF lA 
detector segment. 



^ The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy Observatories, which is operated by 
the Association of Universities for Research in Astronomy, Inc. , (AURA) under cooperative agreement with the National Science Foundation. 
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Because the FUSE SiC channels he on the sun- 
illuniinated side of the spacecraft, their spectra may be 
contaminated by solar emission (Shelton et al. 2001). We 
thus use only data obtained during orbital night to set lim- 
its on emission features in these spectra. In practice, only 
the S VI A933.38 line in our Coma spectrum is affected, as 
all other features in the spectrum are redshifted onto the 
LiF detector channels. All of the Virgo data were obtained 
during orbital night. 

In our Virgo spectrum, the C II* A1037.G2 and O VI 
A1037.62 lines are redshifted to 1040.75 and 1041.36 A, 
respectively, and, if present, might be contaminated by 
flux from the O I airglow features at 1040.94 and 1041.688 
A. To set limits on the flux of C II* and O VI, we include 
the airglow lines in our model spectrum, but limit tlic;ir 
fluxes to that observed in the O I A1027.43 line, which 
is generally much brighter than either of the two longer- 
wavelength features. The resulting limits are higher than 
those for the nearby Ly /3 and O VI A1031.93 lines. Fortu- 
nately, all of the Virgo data were obtained during orbital 
night, when the O I line flux is minimal. 

Finally, we set upper limits on the FUV continuum 
emission from each cluster of galaxies. We consider only 
the region between 1060 and 1070 A, which is free of 
airglow features and for which the LiF lA detector has 
Agff > 20 cm^. For each observation, we determine the 
mean background-event rate from unilluminated regions 
of the detector directly above and below the LWRS ex- 
traction window. Scaling these rates to the size of the 
LWRS window, we find that the observed continuum ex- 
ceeds the detector background by 110 ± 250 and 200 ± 400 
photons cm~^ A~^ sr~^, respectively, for Coma and 
Virgo. (We ignore possible structure in the detector back- 
ground, thought to be present at the 10% level.) 

4. DISCUSSION 
4.1. The Coma Cluster 

From simultaneous fits to EUVE Deep Survey (DS) and 
ROSATV^VC observations of the Coma cluster. Lieu et al. 
(1996a) derive the temperature and emission integral as 
a function of radius for a three-component ICM assum- 
ing a distance of 139 Mpc and an abundance Z = 0.21 
Zq (Hughes et al. 1993). They find a temperature of 
8.7 X 10^ K and an emission integral of 4.8 x 10^^ cm~^ 
for the coolest component within 3' of the cluster cen- 
ter. Combining these parameters with line emissivities es- 
timated with the CHIANTI software package (Dere et al. 
1997) using the Feldman et al. (1992) solar abundances 
scaled to the assumed cluster metallicity, we derive an ex- 
pected O VI A1032 flux of 3.8 x 10"^ photons cm^^ ^-i 
If the O VI emission were evenly distributed within this 
region (an assumption consistent with both EUV and soft 
X-ray maps of the cluster; Bonamente 2000; Briel et al. 
1992), then its surface brightness would be 1600 photons 
cm^^ s^-'^ sr^-'^, slightly less than our 3-cr upper limit to 
the dereddened O VI surface brightness. 

The O VI line emissivity depends critically on the gas 
temperature. From our limit to the O VI flux of Coma, we 
derive an upper limit to the emission integral as a func- 
tion of temperature. In Fig. 2, we compare this limit with 
the best-flt value of the emission integral derived by Lieu 
et al. (1996a) for the central region (r < 3') of the cluster, 



scaled to the area of the LWRS slit. This value of the emis- 
sion integral, EI — 4.2 x 10^^ cm""^, lies just below the 
FUSE 3-a upper limit. If we instead use a 2-cr limit to the 
O VI flux, we can nominally exclude the warm-gas model 
of Lieu et al. (1996a). More recent models of the Coma 
EUV excess do not require a sub- 10^ K gas. Bonamente 

(2000) fits the data with a two-component ICM, the cooler 
of which has a temperature of 2.9 x 10^ K. We would not 
expect significant O VI emission from such a hot gas. 

4.2. The Virgo Cluster 

Two groups (Lieu et al. 1996b; Bonamente et al. 2001) 
have modeled the EUV excess in the Virgo cluster with a 
warm component of the ICM. Both find that its tempera- 
ture rises with radius, while its emission integral falls. The 
brightest FUV emission would thus be expected from the 
cluster center. Unfortunately, the center of the Virgo clus- 
ter is occupied by the galaxy M87, and its strong stellar 
continuum complicates the FUV spectrum of this region 
(cf., Dixon et al. 1996). We thus consider data taken at 
larger radii, about A'A (2242 s) and 6.'0 (8688 s) from the 
cluster center. Both Lieu et al. (1996b) and Bonamente 
et al. (2001) model the warm ICM using a series of concen- 
tric annuli. Our inner data point samples the 3-5 arcmin 
annulus, while our outer point falls in the 5-7 arcmin re- 
gion. Because 80% of our data comes from the outer annu- 
lus, we compare our observations with model predictions 
for the 5-7 arcmin annulus. We estimate that errors in 
our derived limits to the emission integral resulting from 
this simplification are on the order of the uncertainties in 
the FUSE flux cahbration. 

Our upper limit to the emission integral as a function 
of temperature is plotted in Fig. 3. Also shown are the 
best-fit values of the temperature and emission integral 
of the warm ICM, derived from simultaneous fits to the 
DS and PSPC data from the 5 -7 arcmin annulus of the 
Virgo cluster by Lieu et al. (1996b) and Bonamente et al. 

(2001) , respectively, scaled to the area of the LWRS slit. 
The Lieu et al. model parameters are T = 5.8 x 10"'' 
K, EI = 9.35 X 10^1 cm-3, and Z = 0.454 Zq, val- 
ues firmly excluded by our limit to the O VI fiux. The 
Bonamente ct al. model employs a much warmer gas, with 
T = 9.2 X 10'^ K, EI = 4.72 X 10^'^ cm-^, and Z = 0.5 
Zq; these parameters fall just below the FUSE 3-a upper 
limit. If we instead use a 2-cr limit to the O VI fiux, we 
can nominally exclude this model as well. 

The redshift of M87 (z = 0.00436; Smith et al. 2000) 
is slightly greater than that of the Virgo cluster mean. 
Near the center of the cluster, the X-ray emission is well 
centered on M87, and certain asymmetrical features seem 
correlated with radio emission from the nucleus (Belsole 
ct al. 2001). If, in this region, the cooling flow emits at 
the redshift M87, the O VI A1032 line would be strongly 
attenuated by absorption from Galactic C II A1036.3. An 
upper limit to the emission integral derived from the O VI 
A1038 line would be approximately twice our quoted value. 

5. CONCLUSIONS 

We have obtained FUV spectra of two clusters of galax- 
ies with FUSE, the Far Ultraviolet Spectroscopic Explorer. 
Neither spectrum shows significant O VI AA1032, 1038 
emission at the cluster redshift. Our 2-a upper limits on 
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the O VI A1032 flux from Coma and Virgo exclude all pub- 
lished warm-gas models of the EUV excess in these clusters 
and severely constrain the amount of gas at T < 8 x 10^ 
K in either cluster. 

We thank R. Shelton and E. Murphy for discussions 
of FUSE data analysis and D. Buote for thoughtful com- 



ments on the manuscript. We acknowledge the outstand- 
ing efforts of the FUSE P.I. team to make this mission 
successful. This research has made use of NASA's Astro- 
physics Data System and the NASA/IPAC Extragalactic 
Database (NED), operated by JPL for NASA. This work 
is supported by NASA grant NAG 5-5197. 
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Table 1 



3-cT Upper Limits to Far-UV Emission Lines 





Rest Wavelength 


Coma Limits 


Virgo Limits 


Species 


(A) 


Observed 


Dereddened 


Observed 


Dereddened 


S VI 


933.38 


3200^^ 


3700 


4400 


7600 


C III 


977.02 


3900 


4400 


4800 


7700 


NIII 


989.80 


1600 


1800 


4700 


7400 


Nc VI] 


999.20 


1300 


1500 


4500 


7000 


H I (Ly 13) 


1025.72 


1300 


1400 


2000 


3000 


VI 


1031.93 


1600 


1800 


1700 


2600 


C II* 


1037.02 


1200 


1300 


2400 


3600 


VI 


1037.62 


1500 


1700 


2100 


3200 


Ne V] 


1146.1 


2200 


2400 


2900 


4000 



'^To minimize contamination of the SiC 2A channel by solar emission, this upper 
limit is derived only from Coma data taken during orbital night. (All of the Virgo data 
were taken during orbital night.) 

Note. — Units arc photons cm^^ s^^ sr^^. Upper limits are derived assuming an 
intrinsic line width of 40 km (FWHM) at the cluster redshift. Dereddened limits 
assume a CCM extinction curve with E{B — V) = 0.008 for Coma and 0.030 for Virgo 
and Rv = 3.1. 
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Fig. 1. — FUSE spectra of the Coma (top) and Virgo (bottom) clusters, showing the LiF lA spectrum in the region about redshifted O VI 
AA1032, 1038. The Coma spectrum is binned by 8 detector pixels, the Virgo spectrum by 16. The data are presented as histograms and are 
overplotted by models including O VI at our 3-cr upper limit. The observed continuum level is consistent with the dark-count rate determined 
from unilluminated regions of the detector (see text). Wavelengths between 1045 and 1058 A are contaminated by stray light on the LiF lA 
detector segment; the apparent feature at 1054.5 A in the Coma spectrum is thought to be an artifact of the stray-light stripe. 
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5 6 7 8 9 10 

Temperature (10^ K) 



Fig. 2. — Upper limits to the emission integral {EI = J UpHedV) as a function of gas temperature in the Coma cluster. The solid curve 
is derived from the dereddened 3-a upper limit to the O VI A1032 surface brightness, assuming the distance (139 Mpc) and abundance 
(Z = 0.21 Zq) of Lieu et al. 1996a. The dashed line represents a 2-a limit to the same feature. The data point indicates the best-fit value 
of the emission integral and temperature derived for the central region {r < 3') of the cluster by Lieu et al. 1996a, scaled to the area of the 
FUSE LWRS aperture. The FUSE data nominally exclude the warm-gas model at the 2-sigma level. 
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5 6 7 8 9 10 11 12 
Temperature (10^ K) 

Fic;. 3.- Same as Fig. 2, but for Virgo. Our limit to the emission integral assumes the distance (17.2 Mpc) and abundance (Z = 0.5 
Zq) used by Bonamcnte ct al. 2001. The data points indicate the best-fit value of the emission integral and temperature derived for the 
5-7 arcmin annulus by Lieu et al. 1996b (EUV Model 1) and Bonamente et al. 2001 (EUV Model 2), scaled to the area of the FUSE LWRS 
aperture and our assumed distance. Model 1 is completely excluded by our limit, while Model 2 is nominally excluded at the 2-sigma level. 



